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ABSTRACT 
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<.  ■ 

Pulse  mode  sonar  operation  is  analyzed  under  the  assumption  that  the 
scattering  object  T lies  in  the  far  fields  of  both  the  transmitter  and  the 
receiver.  It  is  shown  that,  in  this  approximation,  the  sonar  signal  is  a 
plane  wave  s(x  • 0j  - t)  near  F,  where  0^  is  a unit  vector  directed  from  the 
transmitter  toward  F,  and  similarly  the  echo  is  a plane  wave  e(x  • 0 - t) 
near  the  receiver,  where  0 is  a unit  vector  directed  from  F toward  the 

I trt 

receiver.^  Moreover,  if  F is  stationary  with  respect  to  the  sonar  system 
then  it  is  shown  that 


e(T) 


Re  < 


iTO) 


1 

s((aj)  T.  (o)0  ,(i)0„  )dL)  > 

I 


where  s(aj)  is  the  Fourier  transform  of  s(t)  and  T^(u)0,o)0j)  is  the  scattering 
amplitude  in  the  direction  0 due  to  the  scattering  by  F of  a time-harmonic 
plane  wave  with  frequency  u)  and  propagation  direction  0q . A generalization 
of  this  relation  is  derived  for  moving  scatterers. 
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SECTION  ONE 


INTRODUCTION  - SONAR  ECHO  PREDICTION 

Sonar  systems  for  locating  and  identifying  objects  in  a fluid  environ- 
ment are  used  by  men  and  animals  of  many  species  including  bats,  birds, 
dolphins  and  other  sea  mammals,  and  fish  [1].  Some  species  of  bats  are 
known  to  employ  sophisticated  frequency-modulated  signals  and  to  be  able 
to  discriminate  accurately  among  similar  objects  solely  by  means  of  their 
echoes  [1].  Human  sonar  systems  are  not  yet  this  effective. 

An  active  sonar  system  consists  of  a sound  source  or  "transmitter"  and 
a sound  detector  or  "receiver."  The  transmitter  radiates  an  acoustic  sign  1 
which  is  scattered  by  objects  in  its  environment.  The  resulting  echoes  are 
recorded  by  the  receiver.  In  the  simplest  sonar  systems  the  directions  of 
arrival,  time  delays  and  Doppler  shifts  of  the  echoes  are  used  to  estimate 
the  directions,  ranges  and  speeds  of  the  scattering  objects.  The  goal  of  a 
sophisticated  system  is  to  classify  objects  into  identifiable  classes  by 
means  of  their  echoes.  This  is  an  inverse  scattering  problem.  An  associated 
problem  is  that  of  waveform  design;  that  is,  the  choice  of  a signal  waveform 
which  optimizes  the  echoes  from  prescribed  classes  of  objects.  Before  either 
of  these  problems  can  be  attacked  the  direct  problem  of  echo  prediction  when 
the  signal  and  scattering  object  are  known  must  be  understood.  The  purpose 
of  this  article  is  to  present  a solution  to  this  echo  prediction  problem. 

The  solution  provides  a first  step  toward  the  understanding  of  the  remarkable 
efficacy  of  animal  sonar  systems  which  could  lead  to  the  design  of  more 


effective  man-made  systems. 
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Sonar  systems  are  normally  operated  in  the  pulse  mode  in  which  a 
sequence  of  equally  spaced  short  pulses  is  emitted.  A second  mode  of 
operation  is  the  CW  (continuous  wave)  mode  in  which  a steady  tone  of  fixed 
frequency  is  emitted.  The  construction  of  pulse  mode  echoes  given  below  is 
derived  from  the  theory  of  CW  mode  echoes.  The  analysis  of  each  of  these 
modes  is  based  on  a boundary  value  problem  for  an  acoustic  potential 
function. 

Sonar  echo  prediction  is  analyzed  below  under  the  following  physical 
assumptions . 

• The  sonar  system  (transmitter  and  receiver)  operates  in  a stationary 
homogeneous  unlimited  fluid  medium. 

• The  system  is  stationary  with  respect  to  the  medium. 

• The  sonar  signals  are  generated  by  conservative  force  fields. 

• The  scattering  objects  are  rigid  bodies. 

• The  transmitter  and  receiver  are  in  the  far  field  of  the  scattering 
objects. 

• The  velocities  of  the  scattering  objects  are  less  than  that  of 
sound  and  are  essentially  constant  during  the  Interval  required  for 
the  sonar  pulse  to  sweep  over  the  object. 

• Secondary  echoes  due  to  the  sonar  system  components  are  negligible. 

• Noise  in  the  medium  is  negligible. 

The  case  of  a scattering  object  which  is  stationary  with  respect  to  the 
sonar  system  is  analyzed  first.  The  analysis  is  based  on  the  usual  linear 
theory  of  acoustics  [7].  The  following  notation  is  used,  x - (Xj ,Xj ,Xj) G 
denotes  spatial  coordinates  of  a Cartesian  system  fixed  in  the  medium.  tG  R 
denotes  a time  coordinate,  C R*  denotes  the  domain  exterior  to  the 
scattering  object  and  F = R*  - denotes  the  object.  The  common  frontier  of 
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r and  Q,  which  describes  the  surface  of  the  object,  is  denoted  by  3P  or 
The  equilibrium  state  of  the  medium  is  characterized  by  its  constant  density 
p^y  sound  speed  c^  and  pressure  Pp . It  will  be  assumed  for  notational 
simplicity  that  = 1 and  c^  = 1 since  this  can  be  achieved  by  a suitable 

choice  of  units. 

The  acoustic  field  generated  by  the  transmitter  is  characterized  by  a 
real-valued  acoustic  potential  function 

(1.1)  u = u(t,x),  t € R,  X € () 


which  satisfies  the  inhomogeneous  d'Alembert  equation 


(1.2) 


- Au  = f (t,x)  for  t 6 R,  X € (2 


Here  A is  the  Laplacian  and  f(t,x)  is  a function,  characteristic  of  the 
transmitter,  which  will  be  called  the  source  function.  It  has  the 
structure 

(1.3)  f(t,x)  = 3V(t,x)/9t  + q(t,x) 

where  V(t,x)  is  a potential  for  the  conservative  forces  acting  in  the 
transmitter  and  q(t,x)  is  the  volume  flow  in  the  transmitter  [7,  p.  280]. 

In  what  follows  f(t,x)  is  assumed  to  be  known. 

Pulse  Mode  Scattering.  It  will  suffice  to  analyze  the  scattering  of  a 
single  short  pulse  of  duration  T which  is  emitted  by  a transmitter  that  is 
localized  near  a point  x^ . Hence  the  support  of  f(t,x)  will  be  assumed  to 
satisfy 

(1.4)  supp  f C {(t,x):  tj  < t < tg  + T and  |x  - Xg(  < 6^} 

where  6j  and  t^  are  suitable  constants  and  |x  - Xj]  is  the  distance  between 
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X and  Xjj . The  acoustic  field  generated  by  f(t,x)  is  characterized  by  a 
potential  u(t,x)  which  satisfies  (1.2),  the  Neumann  boundary  condition 
corresponding  to  a rigid  object 


(1.5)  -|^  = Vu  • V = 0 for  t € R,  X 6 

where  V = v(x)  is  a normal  vector  to  9^2,  and  the  initial  condition 


(1.6) 


u(t,x)  = 0 for  t < tj,  X e 


The  primary  field  of  the  transmitter  is  the  field  Uj(t,x)  generated  by 
f(t,x)  when  no  scattering  object  is  present.  It  is  given  by  the  retarded 
potential 


(1.7) 


f(t- 

x-x' 

»x') 

4ir 

x-x' 

dx'  for  t e R,  X e R^ 


where  dx'  = dxj  dx^  dxj . The  variation  of  Ujj(t,x)  near  the  scatterer  F may 
be  simplified  by  using  the  assumption  that  the  transmitter  lies  in  the  far 
field  of  r.  It  will  be  convenient  to  assume  that  the  origin  of  coordinates 
lies  in  F and  to  let  6 > 0 be  the  smallest  constant  such  that 


(1.8)  F c {x:  |x|  < 6} 

Then  the  far  field  assumption  takes  the  form 

(1.9)  |xo  1 » (So  + <5 
If  is  the  unit  vector  defined  by 


(1.10) 


5 


and  if  [x  - x' | is  developed  in  inverse  powers  of  |X(j|,  it  is  found  that 
lx  - x’ 1 = |x„  + (x'  - Xg  - x) I 

= (|xj^  + 2x„  • (x'  - Xg  - x)  + lx'  - Xq  - xl^)!/^ 

(1.11)  = IxJ  (l-26(,*  (x' - Xg  - x)/lxj  + lx' - x^ -xlVlxJ2)i/2 

= 1^0  I (1-00  * 1^0  I + 0(l/lxol^)) 

= |xol  + 00  • ^ ~ ®0  * (^(1/1xq1),  IxJ  ->  oo 

and  the  error  term  is  uniformly  small  for  all  lx'  - x^  1 < <5o  1^1 

Hence,  (1.11)  may  be  substituted  in  (1.7)  and  gives 

(1.12)  uo(t.x)  = s(x-6|-t|lxol)  +Q(i/[xJ^),  lxJ->« 
uniformly  for  t 6 R,  1x1  < 6 where 


(1.13) 


s(t) 


Ait 


f (0p  • (x'  - Xq)  - T,x')dx'  for  T € R 


•'  Ix'-xo  I560 

If  the  error  term  in  (1.12)  is  dropped  the  primary  field  becomes  a plane 
wave  pulse  propagating  in  the  direction  0^.  The  profile  s(t)  will  be  called 
the  signal  waveform.  Note  that  (l.A)  implies  supp  s C [tg-6g-T,  tp  + 6g]. 
The  acoustic  field  produced  when  a plane  wave 


(1.14) 


Ug (t,x)  = s(x  • 0g  - t) , supp  s C [a,b] 


is  scattered  by  F is  the  solution  u(t,x)  of  the  boundary  value  problem 


3^u 

3P- 


- Au  = 0 for  t € R,  X € 12 


i 


I 


(1.15) 
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The  profile  e(T,0,0(j)  will  be  called  the  echo  waveform.  It  depends  on  the 
direction  of  incidence  0^^  of  the  plane  wave  (1.14)  and  the  direction  of 
observation  0.  The  principal  goal  of  this  article  is  to  calculate  the 
relationship  between  the  echo  waveform  and  the  signal  waveform.  The  calcu- 
lation will  be  based  on  the  theory  of 

CW  Mode  Scattering.  CW  mode  fields  are  generated  by  source  functions 
of  the  form 

(1.20)  f(t,x)  = gj  (x)  cos  (i)t  + gj(x)  sin  tot  = Re  {g(x)e 

where  to  > 0 is  a fixed  frequency  and  g = gj  + ig2 • The  corresponding  CW 
mode  field  has  the  same  time-dependence: 

(1.21)  u(t,x)  = Wj  (x)  cos  tot  + Wj  (x)  sin  tot  = Re  {w(x)e 

where  w - Wj  + Iw^ . u(t,x)  must  satisfy  the  d'Alembert  equation  (1.2)  with 
f(t,x)  defined  by  (1.20)  and  the  Neumann  boundary  condition  (1.5).  The 
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initial  condition  (1.6)  is  not  appropriate  for  CW  mode  fields  and  is  replaced 
by  the  Sommerfeld  radiation  condition  which  guarantees  that  u is  a pure 
outgoing  wave  in  the  far  field.  The  corresponding  boundary  value  problem 
for  the  complex-valued  wave  function  w(x)  is 


(1.22) 


(1.23) 


(1.24) 


It  is  assumed  that 


(1.25) 


Aw  + (D^w  = -g(x)  for  X 6 Q 


9w 

9v 


0 for  X 6 9f2 


9w 
9 |xl 


- iojw  = 0(l/lx|^)  for  |x| 


supp  g C {x:  lx  - X 1 < 6g} 


in  agreement  with  (1.4). 

The  primary  CW  mode  field  (x)  is  the  CW  mode  field  generated  by  g(x) 
when  no  scattering  object  is  present.  It  is  given  by 


(1.26)  w„(x)=^ 


} 


i(jj  I x-x ' 


x-x 


g(x')dx'  for  X e 


x'-x„  <6„ 


'o  1-^0 

The  variation  of  (x)  near  F may  be  simplified  by  means  of  the  far  field 
assumption,  as  in  the  case  of  the  pulse  mode  fields.  Substituting  (1.11) 
in  (1.26)  gives 

(1.27)  w„(x)  = + 0(l/lxj^),  lxJ-«> 

uniformly  for  jx|  <6  where 
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(1.28) 


T(<^0o)  = 


e ^ g(x')dx' 


Mx’-Xo1<6„ 


If  the  error  term  in  (1.27)  is  dropped  the  primary  field  becomes  a CW  mode 
plane  wave  of  frequency  o)  propagating  in  the  direction  0^.  It  will  be 
convenient  to  renormalize  the  primary  field  to 


(1.29) 


w„(x,oj0o)  = (271)-^'' e^‘^®o*^ 


The  CW  mode  field  which  is  produced  when  Wq(x,lj0q)  is  scattered  by  F will 
be  denoted  by  w"^(x,W0(j)  . It  is  the  solution 

(1.30)  w''’(x,(jO0n)  = w.(x,(D0.)  + w'^  (x,w0-),  X e 


sc 


of  the  boundary  value  problem 


(1.31) 


(1-32) 


(1.33) 


Aw"*^  + oj^w"^  = 0 for  X E 


9w 

9v 


= 0 for  X E 


9w 


sc 


9 X 


- iwwg^.  = C?(l/|x|^), 


It  is  shown  below  that  in  the  far  field  of  F the  secondary  field 


w^^(x,u)0g)  is  a diverging  spherical  wave 


(1.34) 


Wgc(x,a)0o)  - ^ 


iu)  X 


4TTlxr  X = 1*1^ 


The  coefficient  T^(u)0,u)0j)  is  called  the  scattering  amplitude  of  F . It 
determines  the  amplitude  and  phase  of  the  CW  mode  echo  in  the  direction  0 
due  to  the  primary  wave  (1.29)  with  frequency  w and  propagation  direction  0p . 
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The  solution  of  the  echo  prediction  problem  for  stationary  scatterers 
can  now  be  formulated.  It  is  given  by  the  integral  relation 


(1.35) 


e(T,0,0.)  = Re  { 


e^™  s(o))  T_j^(u)0,w0Q)d(jj} 


where  s(o))  is  the  Fourier  transform  of  the  signal  waveform: 


(1.36) 


S(0)) 


1 

(2T7)  1/2 


-iwT  . V , 
e s(T)dT 


The  derivation  of  this  relation  is  the  principal  result  of  this  article. 

It  is  clear  that  (1.35)  is  not  valid  for  moving  scatterers  because  it 
predicts  that  a signal  and  its  echo  contain  the  same  frequencies,  while  the 
Doppler  shift  in  echoes  from  moving  scatterers  is  well  known.  However, 

(1.35)  can  be  generalized  to  the  case  of  a moving  scatterer  by  first  applying 
it  in  a coordinate  frame  which  is  fixed  with  respect  to  the  scatterer  and 
then  passing  by  means  of  a Lorentz  transformation  to  a frame  which  is  fixed 
with  respect  to  the  sonar  system.  This  calculation  is  carried  out  in  §5 
below. 

Relation  (1.35)  reveals  the  fundamental  importance  for  sonar  echo 
analysis  of  the  scattering  amplitude  T^(o)0 ,0)0^ ) . One  simple  consequence 
of  (1.35)  is  the  result  that  the  echoes  of  high  frequency  pulses  are 
undistorted.  Indeed,  it  lias  recently  been  shown  by  A.  Majda  [3]  that 

(1.37)  11m  T^(oj0,io0o)  = t"(0,0o) 

exists  when  F is  a smooth  convex  object.  It  follows  from  (1.35)  that  if 
s((d)  is  concentrated  in  a high  frequency  band  u)  > Wq  where  T_j^((i)0 ,(i)0p ) is 
nearly  constant  then 


i4i 
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(1.38)  ed.e.e^)  ~ T^(0,e„)Re{  S(u))du)}  = | i;^(e , 0„)  s(t) 

'a  1 / 


by  Fourier's  theorem.  More  interesting,  however  is  the  possibility  of 
applying  (1.35)  to  the  study  of  pulse  distortion  by  reflection  and  its 
dependence  on  the  geometry  of  the  scatterer.  In  particular,  (1.35)  shows 
that  resonances  in  T_^(u)0 ,(i)0(j ) will  produce  selective  enhancement  of  POuaC 
echoes  at  the  resonance  frequencies.  Another  possible  application  of  (1.35) 
is  to  use  Fourier  analysis  of  observed  sonar  echoes  produced  by  known 
signals  to  estimate  the  function  T^(w0,a)0p)  over  a range  of  frequencies. 
These  studies  are  clearly  applicable  to  the  analysis  and  design  of  sonar 
systems . 

This  article  is  based  on  the  author's  monograph  [10]  to  which  reference 
is  made  for  analytical  details.  The  theory  of  scattering  developed  in  [10] 
is  applicable  to  a large  class  of  objects  F with  Irregular,  non-smooth 
boundaries.  The  class  is  defined  by  a "local  compactness  property"  of  the 
domains  (2  = - F [10,  p.  60].  A rather  general  geometric  criterion  for 

the  validity  of  this  property  is  given  in  [10,  p.  63].  The  class  of 
allowable  objects  F includes  all  the  simple  objects  that  arise  in  applica- 
tions such  as  polyhedra,  finite  sections  of  cylinders,  cones,  spheres,  disks 
and,  more  generally,  all  objects  with  piecewise  smooth  surfaces  which  have 
a finite  number  of  smooth  edges  and  vertices.  It  is  assumed  throughout  this 
article  that  has  the  local  compactness  property  so  that  the  results  of 
[10]  are  applicable. 

The  remainder  of  the  article  is  organized  as  follows.  §2  reviews 
briefly  the  definition  and  structure  of  CW  mode  wave  fields.  §3  reviews  tlie 
definition  of  pulse  mode  sonar  wave  fields  and  their  representation  by  means 
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of  CW  mode  fields.  §4  describes  the  theory  of  asymptotic  wave  fields,  as 
developed  In  [10],  and  applies  It  to  the  derivation  of  the  far  field  form 
(1.19)  of  sonar  echoes  and  the  Integral  relation  (1.35).  §5  contains  the 
generalization  of  (1.35)  to  the  case  of  moving  scatterers. 
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SECTION  TWO 

THE  STRUCTURE  OF  CW  MODE  SONAR  ECHOES 

The  definition  and  basic  properties  of  the  CW  mode  fields  w^(x,p)  are 
reviewed  in  this  section.  The  definition  is  based  on  the  boundary  value 
problem  (1.30)  - (1.33)  in  a formulation  appropriate  for  domains  with  non- 
smooth boundaries.  In  the  classical  theory  due  to  Kupradze  [2]  and  Weyl  [8] 
a solution  is  a function  in  the  class  C^(fi)  ri  C*((i)  where  SI  = U u It 

is  known  that  if  9S2  is  smooch  then  a unique  classical  solution  w^(x,p) 
exists  for  each  p = U)6g  S R^ . However,  if  9S^  has  edges  and/or  vertices  then 
at  such  points  the  normal  vector  V is  undefined  and  the  boundary  condition 
(1.32)  is  not  meaningful.  In  this  case  it  was  discovered  that  if  (1.32)  is 
enforced  only  at  points  where  9S2  is  smooth  then  all  solutions  have  first 
derivatives  which  are  singular  at  the  edges  and  vertices  and  no  classical 
solution  exists  [4,  5].  Moreover,  if  the  strength  of  these  singularities 
is  not  controlled  then  the  uniqueness  of  the  solution  is  lost  [5]. 
Physically,  non-uniqueness  can  occur  because  line  or  point  sources  can  be 
situated  in  the  edges  or  vertices.  This  type  of  non-uniqueness  is 
eliminated  by  the  "edge  condition"  [5]  which  requires  that  the  energy  in 
each  bounded  portion  of  (2  should  be  finite.  It  may  be  formulated  mathemat- 
ically as  the  condition  that  u(x)  = w^(x,a)0j)  should  satisfy  [5] 

(2.1)  I { I Vu(x)  I ^ + 0)^  |u(x)  I ^}dx  < °°  for  every  cube  K C R^ 

Kr'ifi 


?RECat)I^O 


page  ELANK-NCT  fILMED 
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In  the  theory  developed  In  [10]  the  edge  condition  (2.1)  is  formulated 
by  means  of  the  function  classes 

(2.2)  L^^((2)  = {u:  u € L2(K  n fi)  for  every  cube  K c R^} 

(2.3)  = L2°‘^(5T)  n {u:  L^‘"(fi)  for  j = 1,2,3} 

j 

where  L2 (M)  denotes  the  usual  Lebesgue  class  of  a set  M C . Condition 
(2.1)  is  equivalent  to  the  condition  uG  L2’^°^(ft).  Moreover,  if 
u = e L2*^°^(^^)  i'  a weak  solution  of  (1.31),  in  the  sense  of  the  theory 
of  distributions,  then  u is  also  in  the  class 

(2.4)  L^^°‘^(A,iT)  = L2**'°‘'(jT)  n {u:  Au  € L2°‘^(i^)} 

For  general  domains  OCR*  and  functions  u € L2*^°^(A,0)  the  Neumann 
condition  (1.32)  is  replaced  by  the  generalized  Neumann  condition  of  [10]: 

(2.5)  I {(Au)v  + Vu  • Vv}dx  = 0 

JO 

for  all  V € L2(0)  such  that  Vv  € L2 (0)  and  v(x)  = 0 outside  of  a bounded 
set.  The  terminology  is  justified  by  the  fact  that  (2.5)  implies  that  the 
classical  Neumann  condition  9u(x)/9v  = 0 is  satisfied  at  every  boundary 
point  X where  90  is  smooth  [10,  p.  41]. 

The  CW  mode  field  w"*"(x,p)  corresponding  to  a primary  field 

(2.6)  W|j(x,p)  = (2it)“^^^  e^*’  *,  p e R* 

and  an  exterior  domain  OCR*  is  defined  in  [10]  by  the  conditions 

(2.7)  w‘^(*,p)  6 L^’^°‘^(A,0)  - L*’*'°‘^(A,0)  n (u:  u satisfies  (2.5)} 
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(2.8) 


3w 


(2.9) 


sc 


9 X 


(A  + |p1^)  w^(x,p)  = 0 for  X G (2 

ilp|  Wgc  ' I^I  “«  where  = w"^  - 


In  this  formulation  the  edge  condition  and  the  generalized  Neumann  condition 
are  incorporated  into  condition  (2.7).  Moreover,  since  it  is  known  that 
weak  solutions  of  (2.8)  are  necessarily  analytic  functions  of  x e (2  [6], 


(2.9)  is  meaningful  for  functions  which  satisfy  (2.7)  and  (2.8).  The 

. +/ 


uniqueness  of  the  functions  w (x,p)  was  proved  in  [10]  for  arbitrary  exterior 
domains  and  their  existence  was  proved  for  domains  f2  which  have  the  local 
compactness  property. 


The  remainder  of  this  section  presents  a calculation  of  the  far  field 


form  (1.34)  of  the  CW  echo  w^^(x,p).  It  is  based  on  the  integral 


representation  [10,  p.  90] 


lx' l=ro 


(2.10) 


X > r. 


where  r^  > 6 and 


(2.11) 


G»(x-x',p)  = 


,1|p1 Ix-x' I 


4tt  fx-x' 


The  as3^ptotlc  form  of  w^^(x,p)  for  x = |x|0,  jx|  may  be  calculated  from 


(2.10)  by  using  the  approximation  |x-x'l  = jx]  - x'  • 0 + 0(1/ |x|)  as  in  §1 . 


1 


The  result  is 


(2.12) 


i|p| |x| 


T^(  |p|9.p)  + I’ 


uniformly  for  all  directions  6 and  p in  any  bounded  set,  where  the  cross 
section  T_j_(p,p')  may  be  written 

f r +/  IN  9 -ip*x  -ip*x  9w^(x,p')i 
(2.13)  T_j^(p,p  ) = J iw  (x.p')  e - e ^ — } dS 


It  follows  from  Green's  theorem  that  the  integral  in  (2.13)  is  independent 
of  rj  provided  that  r^  > 6 and  [p' ] = Ip]. 

The  well  known  reciprocity  theorem  for  T^(p,p')  states  that 


(2.14) 


Vp,p')  = T^(-p'.-p) 


for  all  p and  p'  such  that  jp|  = |p'l.  This  relation  may  be  verified  in  the 
present  setting  by  applying  Green's  theorem  to  w^(x,p)  and  w^ (x,p')  in  the 
region  firj  = (2  O {x:  jxj  < ro)  and  using  the  representation  (2.13).  If 

is  not  smooth  the  calculation  must  be  based  on  the  version  of  Green's 
theorem  given  by  the  generalized  Neumann  condition  (2.5).  The  technique 
for  doing  this  is  described  in  [10,  pp.  57-8]. 

It  can  be  shown  that  if  the  scatterer  F has  a piecewise  smooth  surface 
9F  = 9fi  then  w"*" (x,p)  is  continuous  and  9w'*'(x, p) /9v  = 0 at  all  boundary 
points  except  edges  and  vertices.  In  this  case  Green's  theorem  can  be 
applied  to  Wo(x,-p)  and  w"*"(x,p')  in  (2(rj)  . This  calculation,  together  with 
(2.13),  implies  the  representation 


T+(p,p')  = I w'‘’(x,p')  ^ (e  ’') 


Moreover,  under  the  same  hypothesis  the  application  of  Green's  theorem  to 


Gq(x-x',p)  and  w^(x',p)  gives  the  representation 


(2.16)  w''’(x,p)  = w (x,p)  + [ w~''(x' ,p)dS ' , x6 

ha 

and  the  unknown  boundary  values  of  w"*"(x,p)  satisfy  the  Fredholm  equation  [4] 

(2.17)  w‘*’(x,p)  = 2w.(x,p)  + 2 w'^(x',p)dS'.  x € 

Jail 

These  relations  provide  a basis  for  both  exact  and  approximate  calculations 
of  w"'’(x,p)  and  T (p,p'). 
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SECTION  THREE 

THE  STRUCTURE  OF  PULSE  MODE  SONAR  ECHOES 

In  this  section  the  Hilbert  space  method  of  [10]  is  used  to  construct 
the  sonar  echoes  produced  by  pulse  mode  signals.  The  construction  is  based 
on  the  fact  that  the  negative  Laplacian  A = -A,  acting  on  functions  that 
satisfy  the  edge  condition  and  generalized  Neumann  condition,  defines  a 
selfadjoint  operator  in  the  Hilbert  space  L2(£^).  It  was  shown  in  [10]  that 
the  CW  mode  fields  w^(x,p)  form  a complete  set  of  generalized  eigenfunctions 
for  A.  This  result  is  used  below  to  derive  an  integral  representation  of 
pulse  mode  sonar  echoes. 

The  edge  condition  and  generalized  Neumann  condition  for  function  in 
are  formulated  in  [10]  by  means  of  the  function  classes 

(3.1)  L*(fi)  = LjCS))  n {u:  LjCfi)  for  j = 1,2,3} 

(3.2)  L*(A,fi)  = L^(S))  n {u:  Au  e 

(3.3)  L2(A,fi)  = L2(A,fi)  n {u:  u satisfies  (2.5)  for  all  v e L2(fi)} 

The  operator  A in  L2(f2)  is  defined  by 

(3.4)  D(A)  = L2(A,0) 


(3.5)  Au  = -Au  for  all  u € D(A) 

It  was  shown  in  [10,  p.  41]  that  A is  a selfadjoint  non-negative  operator 
In  L (U): 
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(3.6) 

A = A*  > 0 

The  non-negative 
It  was  shown  in 

square  root  of  A,  denoted  by  A*^^,  is  also  needed  below 
[10]  that 

(3.7) 

D(a‘/")  = L^(fi) 

(3.8) 

llA*/^u«  = IIVull  for  all  u E D(A*/^) 

where  II  *11  denotes  the  L2(fl)-norin. 

The  d'Alembert  equation  (1.2)  is  interpreted  below  as  an  ordinary 
differential  equation 

(3.9)  + Au  = f(t,‘),  t e R 

for  a function  t -*■  u(t,*)  e L2(fJ).  It  is  assumed  that 

(3.10)  supp  f C {t:  to  < t < tjj  + T} 

and  a solution  of  (3.9)  is  sought  which  satisfies 

(3.11)  u(t,*)  = 0 for  t < to 
in  agreement  with  (1.4),  (1.6). 

A formal  solution  of  (3.9),  (3.11)  is  defined  by  the  Duhamel  integral 


(3.12) 


u(t,-) 


ft 

sin  (t  - T)A^^^}f  (t,  *)dT,  t > t 
J t “ 


where  A"'^^  sin  t A^^^  is  defined  by  means  of  the  spectral  theorem.  Integra- 
tion by  parts  in  (3.12)  gives  the  alternative  representation 


i 

't 

\ 


% 
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u(t,*)  = I {cos  (t-T)A'^^}  j f(t',0  dt'  dT 

In  tf) 


It  will  be  assumed  that  the  source  function  f satisfies 


(3.14) 


f e C(R,L*(J1))  = C(R,D(A^/^)) 


It  can  then  be  verified  that  (3.12)  defines  a function 


(3.15) 


u6  C^(R,L,(i2))  O C^(R,D(A‘(^))  n C(R,D(A)) 


whose  first  and  second  t-derivatives  satisfy 


(3.16) 


^ = j {cos  (t-  r)A^^^}  f(T,*)dT 


(3.17) 


= f(t 


.•)  - P {A^/^ 


sin  (t  - t)A*/^}  f (t,  •)  dT 


= f (t, •)  - Au 


In  particular,  u(t,*)  € D(A)  for  every  t e R and  hence  u(t,x)  satisfies  the 
edge  condition  and  generalized  Neumann  condition  for  every  t.  The  uniqueness 
of  this  kind  of  solution  was  proved  in  [9]. 

The  remainder  of  this  article  deals  with  acoustic  fields  of  the  form 
(3.12)  at  times  t > t^  + T;  l.e.,  after  the  primary  field  has  been 
established  and  the  sources  have  ceased  to  act.  It  will  be  assumed  that 


(3.18) 


f(t,x)  dt  “ 0 for  all  x € 
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This  hypothesis  is  made  to  simplify  the  analysis  and  is  not  essential  to 
the  method.  The  physical  meaning  of  (3.18)  is  clear  from  (1.3). 
Equations  (3.13)  and  (3.18)  imply  that 


u(t,*)  = 


(3.19) 


to+T 


{cos  (t  - t)A*^^} 


f(t',*)dt'dT 


for  t > tg  + T 


Moreover,  since  f(t,x)  is  real-valued  (3.19)  implies  that 


(3.20) 

where 


u(t,x)  = Re  {v(t,x)}  for  t > t^^  + T 


(3.21) 


v(t,-)  = 


and 


(3.22) 


h = 


to+T 


iTA^/2 

e 


f(t',*)  dt'dT 


These  equations  contain,  in  abstract  form,  the  solution  of  the  sonar  echo 
prediction  problem.  A more  concrete  representation  is  provided  by  the 
eigenfunction  expansion  theorem  of  [10]  which  will  be  reviewed  here  briefly 
and  applied  to  (3.21). 

The  expansion  theorem  states  that  every  h e Lj  (fi)  has  a generalized 
Fourier  transform 


(3.23)  h.(p)  “ L2(R^)-llni  f w"^(x,p)  h(x)  dx 

M-~" 

where  r'»  (x:  |x|  < m)  and  a corresponding  eigenfunction  representation 
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for  every  bounded  measurable  function  4'(X)  which  is  defined  for  X > 0.  In 
what  follows  the  relations  (3.23),  (3.24)  will  be  written  in  the  symbolic 
form 

(3.27)  h,(p)  = w"'’(x,p)  h(x)  dx 

(3.28)  h(x)  = I w"'’(x,p)  h^(p)  dp 

However,  theae  integrals  are  not  in  general  convergent  and  must  be 
Interpreted  in  the  sense  of  (3.23),  (3.24). 

A second  complete  family  of  generalized  eigenfunctions  for  A is  defined 
by 

(3.29)  w~(x,p)  = w"'’(x,-p) 


This  family  can  also  be  characterized  by  the  boundary  value  problem  (2.6)- 
(2.9),  but  with  the  outgoing  radiation  condition  (2.9)  replaced  by  the 
incoming  radiation  condition 
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The  relations  (3.23)  - (3.26)  hold  with  w (x,p),  h^(p)  and  replaced  by 

/N 

w~(x,p),  h (p)  and  respectively.  This  was  proved  in  [10].  It  also 
follows  directly  from  (3.23)  - (3.26)  and  (3.29).  The  far  field  behavior 
of  w (x,p)  follows  from  (2.12)  and  (3.29): 


(3.31) 


Wsc(x,p) 


g-i|pl  Ul 

4tt1x| 


T_(  Ip  l0,p)  + 


00 


where 


(3.32)  T_(p,p')  = T^(p,-p') 
Application  of  (3.26)  to  (3.21)  gives 

(3.33)  (<I>^  e“^*^  ^ h)(p)  = h^(p) 

and  hence,  by  (3.24) 


(3.34) 


v(t,x) 


w'*'(x,p)  e h_^(p)  dp 


The  family  {w~(x,p):  p £ R^}  provides  a second  representation 


(3.35) 


where 


(3.36) 


v(t,x) 


w”(x,p) 


h_(p)  dp 


R 


3 


h_(p) 


w (x,p)  h(x)dx 


If  h is  defined  by  (3.22)  where  f satisfies  (3.10),  (3.14),  (3.18)  then  the 
continuity  of  and  together  with  (3.26),  imply  that 
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^ 1 <Mt  ± 

(3.40)  f+(w,p)  = . 1/2  e w (x,p)  f(t,x)dxdt 

■'to 

The  total  field  produced  by  the  scattering  of  acoustic  waves  with  a 

/N 

source  function  f(t,x)  is  represented  by  (3.34)  and  (3.35)  where  h^(p)  and 
h_(p)  are  defined  by  (3.39),  (3.40).  These  equations  are  not  directly 
applicable  to  the  boundary  value  problem  (1.14)  - (1.17)  for  the  scattering 
of  plane  waves.  However,  it  will  be  shown  that  this  problem  can  be  reduced 
to  the  preceding  one.  To  this  end  let  to  and  R be  radii  such  that  6 < r^  < R 
and  let  j (x)  G C (R  ) be  a function  with  the  properties 

(3.41)  j(x)  = 0 for  |x|  < rg,  j(x)  = 1 for  |x|  > R 

The  solution  of  the  boundary  value  problem  (1.14)  - (1.17)  will  be  constructed 
as  a sum 


0.42) 


u(t,x)  = j(x)  Ug(t,x)  + u^^(t,x) 


Note  that  u^^(t,x)  coincides  with  the  echo  (1.18)  for  |x  | > R and  all  t € R, 
by  (3.41)  . 

The  solution  of  (1.14)  - (1.17)  is  given  by  (3.42)  if  u^^(t,x)  satisfies 
9^u 

(3.43)  ^ - Au  = f(t,x)  for  t € R,  x € 

3t^ 

3u 

(3.44)  — — = 0 for  t e R,  X e 

3v 

(3.45)  u (t,x)  = 0 for  t+b+R  <0,  x € 

sc 

where 

(3.46)  f(t,x)  = - - A j(x)  Ug(t,x) 

Note  that 

(3.47)  supp  f C {(t,x):  -b  - R < t f -a+R  and  < jx]  5 R} 

by  (1.14)  and  (3.41).  In  particular,  (3.10)  holds  with  t^  = -b  - R and 
T = 2R  + b - a.  Moreover,  (3.14)  holds  if  s(x)  € L2(R)  = L2(R) 
fi  {s(t):  s'(t)  and  s"(t)  are  in  L2(R)},  and  (3.18)  holds  if 

(3.48)  [ s(T)  dx  = 0 

■'a 

Of  course,  if  3F  is  not  smooth  then  u^^  must  satisfy  the  edge  condition  and 

~ N 

generalized  Neumann  condition;  i.e.,  u (t,*)  € D(A)  = L (A,fi)  for  all  t G K. 

SC  2 

The  solution  is  given  by 

(3.49)  Ug^(t,x)  = Re  (Vg^(t,x)} 


where  v^^  has  the  two  representations 
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(3.50) 


V (t 

sc' 


,x)  = j w^(x,p)  e h^(p)dp 


corresponding  to  w and  w , and  h^(p)  and  h_(p)  are  defined  by  (3.39),  (3.40) 

/s 

and  (3.46).  This  section  is  concluded  with  a calculation  of  h_(p) . 

The  support  of  f in  space-time  is  bounded,  by  (3.47),  and  a simple 
calculation  gives 


(3.51) 


e f(j.^x)dt  = s(-w)  (A+O)^)  j(x)  e * 


where  s(o))  is  defined  by  (1.36).  Thus,  by  (3.29)  and  (3.51), 


(3.52) 


f_(u),p)  = s(-(u) 


f w+(x,-: 

Jo 


p)(A+w^)  j(x)  e ^ dx 


Now  by  (3.41) 


(3.53)  supp  (A+to^)  j(x)  e ^ C {x:  Tq  5 !*l  < R)  = 

Hence,  replacing  by  in  (3.52)  and  applying  Green's  theorem  to 

K rj 

w^(x,-p)  and  j (x)  ^ in  the  region  r^  < |x|  < R gives 


J.HpI.p)  - 3<|p1)  |(«^-.-p)  ■*= 


(3.54) 


= s(1p|)  T^(-|p|0o,-p)  = s(1p|)T^(p, |p|6o) 


by  (2.13)  and  (2.14).  Thus,  by  (3.39) 


(3.55) 


h_(p)  = (2ir)‘/2i  |p|-i  S(|p|)  T^(p,|pl0„) 


Combining  (3.49),  (3.50)  and  (3.55)  gives  the  representation 


1 
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(3.56)  Ug^(t,x)  = Re<|(2TT)‘^^  i | w“(x,p)e  ^ ^ ^ | p | "‘s(  | p | )T^(p,  | p 1 6g ) dp|> 

R* 

for  |x|  > R.  Actually,  (3.56)  is  valid  for  all  x £ 0.  This  can  be  seen  by 
noting  that  (3.56)  does  not  depend  on  j (x)  and  that  condition  (3.41)  may  be 
replaced  by  the  condition  that  supp  V j (x)  is  contained  in  an  arbitrarily 
small  neighborhood  of  F. 


29 


SECTION  FOUR 

SONAR  ECHOES  IN  THE  FAR  FIELD 

The  sonar  echo  u^^(t,x)  of  the  primary  field  Up(t,x)  = s(x  • 6q  - t) 
originates  at  F and  reaches  points  x in  the  far  field,  characterized  by 
jx]  » 5,  after  a time  interval  of  magnitude  comparable  with  jx].  Hence  the 
far  field  form  of  u^^(t,x)  coincides  with  its  asymptotic  form  for  large  t. 

The  latter  is  provided  by  the  theory  of  asymptotic  wave  functions  developed 
in  [10].  The  relevant  parts  of  the  theory  are  summarized  here  and  applied 

to  Ug^(t.x) . 

The  simplest  theorem  concerning  asymptotic  wave  functions  states  that  if 

-itA^/2 

(4.1)  u(t,x)  = Re  {v(t,x)},  v(t,0  = e h,  h e LjCii) 

then 

(4.2)  11m  Du(t,*)  - u°°(t,*)ll  = 0 

t-K» 

CO 

where  u (t,x),  the  asymptotic  wave  function  for  u(t,x),  is  defined  by 

(4.3)  u“(t,x)  = 
and 

(4.4)  F(t,6)  = Re  <[-^—375-  e^™  h_(u)0)  (-iu))do) 

where  h_(p)  is  defined  by  (3.36).  This  result  Is  proved  in  [10], 


Corollary  8.3. 
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Equations  (4.1)  - (4.4)  can  be  applied  to  the  echo  u^^(t,x) . Indeed, 
(3.56)  Implies  that 

1/2 

(4.5)  Ug^(t,x)  = Re  {Vg^(t,x)},  Vg^(t,*)  = e h 

where  h Is  defined  by  (3.55).  The  calculation  leading  to  (3.55)  shows  that 
h Is  also  given  by 


(4.6) 


rto+T  .T 


f (t,x)dtdT 


where  f(t,x)  Is  defined  by  (3.46).  A simple  calculation  based  on  (3.46) 
and  (4.6)  gives  the  alternative  representation 


(4.7)  h(x) 


e < 


|Aj(x)  j" 


s(t)dt  - 6o‘^J  s(x*  0j-T)^dT 
X*0o-T 


It  follows  that  If 


(4.8) 


s e Lj(R),  supp  s c [a 


.b],  r 

» S3 


s(t)dT  = 0 


as  was  assumed  above,  then  he  D(A) . In  particular,  he  Lj  (S7)  and  u^^(t,x) 
has  the  asymptotic  wave  function 


(4.9) 


e(| 

1x1 

-tj 

,0.0o) 

where  e(T,0,0Q)  = F(t,0)  Is  given  by  (4.4)  with  h_  defined  by  (3.55) 
Substitution  gives  equation  (1.35)  of  §1;  l.e.. 


(4.10) 


e(T,0,0o)  - Re<  e s(o))  T 


_j_(u)0,u)0|,)da)|> 


Moreover,  (1.19)  of  §1  now  has  the  Interpretation 
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(4.11) 


t->oo 


Stronger  convergence  results  follow  from  the  hypotheses  (4.8).  Indeed, 
(4.8)  implies  that  h £ D(A*/^)  = L2(51)  and  hence  Theorem  8.5  of  [10]  implies 
that 


(4.12) 


where  x.  = t. 


(4.13) 


lim 

t“*00 


9usc(t,‘) 


9x. 

J 


u .(t,*) 
sc.j' 


= 0,  j = 0,1, 2, 3 


e.(|x|-t,0,eo) 

u .(t,x)=^ 


, j = 0,1, 2, 3 


e (x,0,eo) 


_ 9e(T,0,0o) 
9t 


(4.14) 


and 


(4.15) 


Note  that 


(4.16) 


= Re  { j e^™  (-ill))  s(o))  T^(aj0,(i)0p)dw} 

* ft 


ej(T,0,0o)  = -0j  ej(T,0,0o)  for  j = 1,2,3 


u (t,x)  = 
sc , 0 ’ 


9u  (t,x) 


sc 


9t 


The  remaining  asymptotic  wave  functions,  corresponding  to  the  space  deriva- 
tives, are  not  the  derivatives  of  u°°  but  differ  from  them  by  terms  which 
converge  to  zero  In  L2(^)  when  t o“.  This  can  be  proved  by  means  of  Lemma 
2.7  of  [10]. 
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SECTION  FIVE 

SONAR  ECHOES  FROM  MOVING  OBJECTS 

The  preceding  sections  treated  scattering  from  objects  that  were 
stationary  with  respect  to  the  sonar  system.  In  this  section  the  analysis 
is  extended  to  moving  scatterers  under  the  simplifying  assumption  that 
their  velocities  are  subsonic  and  essentially  constant  during  the  scattering 
of  a pul se . 

The  relative  velocity  of  the  sonar  system  and  the  scatterer  during  the 
scattering  will  be  described  by  a constant  vector  v with  magnitude  v = |v| 
such  that  0 < V < 1 ( = speed  of  sound).  The  space-time  coordinates  of  a 
Galilean  reference  frame  in  which  the  sonar  system  is  at  rest  will  be 
denoted  by  (t,x)  = (t,Xj ,X2,X3) , as  in  the  preceding  sections.  The  coor- 
dinates of  a second  Galilean  frame  in  which  the  scatterer  F is  at  rest  will 
be  denoted  by  (t',x')=  ( t ' , x| , x^ t xp . The  two  coordinate  systems  are 
related  by  a Lorentz  transformation  based  on  v.  It  may  be  assumed  that  the 
spatial  axes  of  the  two  systems  are  parallel  and  that  the  origin  O'  of  the 
second  system  moves  along  the  Xj-axis  in  the  positive  direction  with  speed 
V.  The  Lorentz  transformation  then  takes  the  form 
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(5.2) 


t'+vx| 

/l-v^ 


x[+vt' 


The  assumption  that  F is  in  the  far  field  of  the  transmitter  implies 
that  the  sonar  signal  is  effectively  a plane  wave  near  F,  as  was  shown  in 
§1.  Thus  the  signal  is  described  relative  to  the  sonar  system  frame  by  a 
scalar  potential 


(5.3)  Uj (t,x)  » s(x  • 00  - t) 

The  Lorentz  transformation  (5.1)  satisfies  the  relation 

(5.4)  X • 0„  - t = Yo(x'  • 0J  - t') 


where 

(5.5) 

and 


1-V0‘  l-v-00 

/l-v^  /1-v^ 


Ii  _ 


1-V0[ 


(5.6)  0 

Thus  the  signal  takes  the  form 


„,2  /1-V^0^  ,,3  _ /l-V^0Q^ 

0 " 1-V0^  ’ 0 1-V0J 


(5.7)  uj(t'.x')  = s’(x'  • 0^  - t’) 

relative  to  the  object  frame  where 


(5.8)  s' (t)  = s(YqT) 

It  is  easy  to  verify  that  (5.6)  defines  a mapping  of  the  unit  vectors  Oq  € r’ 
into  unit  vectors  0J  € R*. 
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The  d'Alembert  equation  (1.15)  is  invariant  under  the  Lorentz  trans- 
formation (5.1).  Hence,  relative  to  the  object  frame  coordinates  (t',x') 
the  echo  produced  by  the  scattering  of  the  signal  (5.7)  from  F is  the  solu- 
tion of  the  same  boundary  value  problem  (1.15)  - (1.18)  that  was  solved  in 
§3.  It  follows  from  the  analysis  of  §4  that  the  echo  is  given  in  the  far 
field  by  the  asymptotic  wave  function 


(5.9) 


u '(t'.x')  = 
sc  ’ 


e'dx' 

-t’,0'.0i) 

[x'l 

X = X 


with  echo  waveform 


(5.10) 


e'(T,e',0')  = Re< 


^iTU)  T_^((jj0  ' ,(i)6j)d(j0 


where  T^(p,p')  is  the  same  function  defined  in  §2. 

Equations  (5.9),  (5.10)  describe  the  echo  relative  to  the  object  frame. 
To  describe  it  relative  to  the  sonar  system  frame  it  is  necessary  to  apply 
the  Lorentz  transformation  again.  This  may  be  done  by  noting  that  along  the 
ray  x'  = | x' [ 0 ' with  0 ' f ixed  and  | x ' | > 0 one  has  j x'  | = x'  • 0 ' . Hence, 
using  (5.4)  - (5.6)  with  0q  replaced  by  0, 


e’dx'l  - 


t’,0',0;) 


(5.11) 


where 


(5.12) 


(5.13) 


e'(x'  • 0'  - t',0',0o') 

e(x  • 0 - t,0,Oo)  = e( |x|  - t,0,Ofl) 


e(x,0,0o)  = e'(Y~  T,0',0o) 

1-V0*  l-v*0 


Y = 


t^l-v^  /l-v^ 


ak 
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and  6j  (respectively  9')  are  derived  from  6,  (respectively  6)  by  (5.6). 
Thus  the  echo  is  described  relative  to  the  sonar  system  frame  by 


(5.14) 


where  the  echo  waveform  e(T,6,0o)  is  defined  by  (5.10),  (5.12)  and  (5.13). 
The  resulting  form  of  e(T,0,0Q)  can  be  written 

)dcoV 


(5.15) 


e(T,9,9o) 


= JL 
Yo 


Re  < 


iTti) 


JL 

,Yo 


T_j_(a)Y0’,o)Y6o 


If  V = 0 then  Y “ Yo  “ 9'  “ Qq  = 9o  (5.15)  reduces  to  equation 

(4.10)  for  the  echo  waveform  of  a stationary  scatterer.  Comparison  of 
(4.10)  and  (5.15)  shows  that  the  motion  of  the  scatterer  relative  to  the 
sonar  system  causes  three  kinds  of  distortion  in  the  echo.  These  are  a 
frequency  shift  in  the  signal  waveform  of  amount 


(5.16) 


Ato  = — 1, 

Yo 


a frequency  shift  in  the  scattering  amplitude  T^  of 


amount 


(5.17) 


A'oj  = Y ~ 1 


and  an  angular  distortion  in  the  scattering  cross  section  T^  due  to  the 
replacement  of  9 and  0^^  by  0'  and  9^.  The  frequency  shift  (5.16)  is  the 
usual  Doppler  shift. 
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